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ABSTRACT 

Studies of a class of infinite one dimensional self-gravitating systems have highlighted 
that, on the one hand, the spatial clustering which develops may have scale invariant 
(fractal) properties, and, on the other, that they display "self-similar" properties in 
their temporal evolution. The relevance of these results to three dimensional cosmo- 
logical simulations has remained unclear. We show here that the measured exponents 
characterizing the scale-invariant non-linear clustering are in excellent agreement with 
those derived from an appropriately generalized "stable-clustering" hypothesis. Fur- 
ther an analysis in terms of "halos" selected with a friend-of-friend algorithm reveals 
that such structures are, statistically, virialized across the range of scales corresponding 
to scale-invariance. Thus the strongly non-linear clustering in these models is accu- 
rately described as a virialized fractal structure, very much in line with the "clustering 
hierarchy" which Peebles originally envisaged qualitatively as associated with stable 
clustering. If transposed to three dimensions these results would imply, notably, that 
cold dark matter halos (or even subhalos) are 1) not well modeled as smooth objects, 
and 2) that the supposed "universality" of their profiles is, like apparent smoothness, 
an artefact of poor numerical resolution. 
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1 INTRODUCTION 

The much acclaimed successes of the ACDM cosmology in 
matching many observations concern essentially its homoge- 
neous limit and the linear regime in which perturbations to 
homogeneity are small. The success of the model in account- 
ing for the numerous observations which probe the non- 
linear regime, where density fluctuations are large, is much 
more uncertain. An important consideration in this respect 
is the great difliculty of calculating the model's predictions 
in this regime. Even in the idealized limit in which clus- 
tering arises from gravity alone, predictions reside solely on 
numerical simulations. The latter have, despite impressive 
increases in their size and sophistication, still quite limited 
spatial resolution. Further, essentially because there are no 
non-trivial analytical benchmarks which they can be tested 
against, it remains unclear to what extent their results are 
conditioned by these resolution limits. In this paper we ex- 
plore what might be learnt about the nature of non-linear 
gravitational clustering from the study of a class of one di- 
mensional (ID) models, which have the interest of offering, 
even with modest computer resources, very much greater 
spatial resolution and exact numerical integration. An early 



study in this spirit is that of iMelottI (|l982l . Il983f ) which 
used such a ID mo del to explo r e clus tering in hot dark 
matter cosmologies. iRouet et al.l (|l990t ) derived and stud- 
ied a slightly different model to the most naive ID version 
of clusteri ng in an Einstein de S i tter (EdS) universe con - 
sidered by fyano fc Goudal (Il998l): lAureU fc Fanelhl (120021): 
lAurell et all (|2003l ). IValagead (|2006l ) and iGabrielli et all 
(|2009^ liave studied the version without expansion, while 
Miller et al. jMiUer fc Rouetll2002l . l2006l : iMiller et al.ll2007l : 
iMiller fc Rouetl 2010l ) have reported results on all three of 
these models. Despite these studies it remains unclear, how- 
ever, whether a good, and really useful, analogy can be 
made wit h 3D cosmolog i cal si mulat ions. Generaliz i ng th e 
resuhs of lYano fc Goudal (|l998l ) and lGabrielh et al.l |200i), 
we confirm here the very strong qualitative similarities in 
the temporal development of clustering to that in the 3D 
case. Further we show the scale invariant properties of 
the clustering in real spac e, emphasized and analyzed by 
iMiUer fc Rouetl (|2002| . |2006| ): IMiller et all (|200'iD . can in fact 
be well explained within an analytical fra mework similar 
to one pro posed long ago for the 3D case (|Peeblesl [l973 : 
iDavis fc Peeblca,1977. : ,Pccblca,1980l ). We derive the appro- 
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priate ID version of this generalized "stable clustering" pre- 
diction for the correlation exponents, and show it to provide 
a very good approximation to the numerical results. We also 
show that, in a precise sense, the strongly non- linear cluster- 
ing can be properly described as a virialized fractal hierar- 
chy, in line with what was originally anticipated qualitatively 
in three dimensions. We argue finally for the extrapolation 
of these conclusions about the qualitative nature of cluster- 
ing to the relevant 3D case, and consider its implications, 
in particular for what concerns the nature and properties of 
"halos" in cold dark matter cosmologies. 



2 ID VERSIONS OF COSMOLOGICAL 
N-BODY SIMULATIONS 

Dissipationles s cosmologic al N-body simulations (for a re- 
view see, e.g. jBaglal |200^) solve numerically the equations 
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where Xi are the comoving particle coordinates, a{t) is the 
appropriate scale factor for the cosmology considered, with 
Hubble constant H{t) — a/a. The superscript 'J' in the sum, 
which runs in practice over the infinite system constituted 
by periodic copies of a cube containing particles, indicates 
that the (badly defined) contribution of the mean mass den- 
sity (which sources the Hubble expansion) is subtracted. We 
consider here simply the ID system obtained by replacing 
the 3D Newtonian force term by the analogous ID expres- 
sion, derived starting from the ID Poisson equation which 
gives a force between particles in one dimension (or, equiva- 
lently, infinite parallel sheets embedded in three dimensions) 
independent of separation. We thus consider the equations 



lim sgn(xi 
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where the limiting procedure used in the sum is simply 
a conve nient way to e xplicit the subtraction of the back- 
ground (|Kiesshng||2003l ). and g is the coupling constant (with 
g = 2-k'EG for sheets of surface mass density E). Initial con- 
ditions are generated, just as in 3D cosmological simulations, 
by applying appropriate small perturbations to particles ini- 
tially on an infinite perfect lattice. If the displacement from 
its initia l lattice site of parti cle i is Ui it can be shown rig- 
orously (|Gabrielli et al.ll2009l ) that the net force on the par- 
ticle, until it crosses another one, is exactly proportional to 
Ui — (u), where (u) is the average displacement. The equa- 
tions of motion up to the time at which particles cross be- 
come 



ili + 2Hui 
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(3) 



where no is the mean particle density (and we take (u) — 0). 
When the particles in our model cross one another the force 
changes simply by ±2g. However, since in ID a crossing of 
particles is equivalent, up to exchange of labels, to an elastic 
"collision" between them, one can instead consider (as we 
are not interested in the labels of the particles) a system 
of particles which bounce elastically and follow Eq Q at 
all times between "collisions" . Finally, by a transformation 



of the time coordinate to r = J dt/a'^^^, and appropriate 
choice of units of r, we can rewrite these equations as 



dr^ dr ' 



(4) 



i.e., as those of a set of damped inverted harmonic oscil- 
lators. We note that these equations coincide simply with 
those for displacements of fluid e l emen ts in the Zeldovich 
approximation (see, e.g. iBuchertI (|l992l )). which is in fact 
exact up to "shell crossing" for ID perturbations. The same 
system can thus be derived as a simple anal ytical continua- 
tion of this approximati on at shell crossing (jYano fc Goudal 
1 19981 : lAurell et al.ll2003l ) . 

For a generic cosmology F in Eq. Q is a non-trivial 
function of t, but in the specific case of an EdS cosmol- 
ogy (for which a oc t^^^) it is a constant, F = l/\/6 (and 
r cc I n a). The model derived and studied bv iRouet et al.l 
(|l990l ) corresponds instead to the case F — 1/V2, while the 
case F = (i.e. static limit) has also been studied by several 
authors (see references given above). We note that, despite 
the fact that we obtained the value F = I/VG above, this 
is not necessarily the more appropriate value to consider for 
our study: in the derivation just given the 3D Hubble law 
has been imposed by hand. This has been done because the 
ID Hubble law (which would arise if one follows fully the 
analogy with the 3D derivation starting from physical co- 
ordinates) gives a completely different qualitative behavior 
— collapse in a finite time — which has no relevance to 3D 
cosmology. Thus in order to "mimic" EdS cosmology, one 
could equally take the functional form of the corresponding 
Hubble law, but leave its normalization as a free parameter. 
This leads to the Eqs. Q with F an undetermined constant. 
These equations then simply define a simple class of toy 
models, which may be useful for qualitative understanding 
of the 3D problem. 

Given initial displacements and velocities, Eqs. Q can 
trivially be integrated exactly between particle collisions 
when F is a constant. To determine the time of the next 
collision (and which pair of particles it involves) requires 
then only the solution of algebraic equations. The numeri- 
cal integration can therefore be perfor med "exactly", i .e., up 
to machine precision. As described bv lNoullez et all (|2003l ) 
the integration can be sped up optimally using a "heap" 
structure. All the results presented here are for systems with 
N — 10'' particles (and periodic boundary conditions). 

In making the analogy with 3D simulations, which aim 
to reproduce the collisionless limit of gravitational cluster- 
ing and include a smoothing of the singularity in the 3D 
Newtonian force for this reason, it is to be noted that, on 
the time scales we will consider, these ID systems, without 
force smoothing, can be expected to represent extremely well 
this limit. Studies of col lisional re laxation in finite I D self - 
gravitating systems (see ljovce fc WorrakitpoonponI (|201(j ). 
and references therein) have shown that it is very sup- 
pressed compared to that in 3D systems, a typical relax- 
ation time for a virialized system of number density n being 
~ (10* — 10^)N/^^gn. Using this estimate in the simulations 
reported below, it is simple to show that even the relaxation 
time for the smallest and densest clusters is much longer 
than the duration of the simulation. The reason for this rel- 
ative suppression of this relaxation is that, in one dimension, 
there is no analogy to 3D 2-body relaxation: the "collisions" 
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of particles we have discussed above do not contribute to 
the usual collision term in the Boltzmann equation. Indeed 
it is for this reason that particle crossings and "collisions" 
are equivalent. 



3 COMPARISON WITH SIMULATIONS IN 
3D: SELF-SIMILARITY 

Various studies (see references above) of these models start- 
ing from close to uniform initial conditions have noted the 
striking qualitative similarity of the clustering to that in 
three dimensions; one observes the formation of structures 
(i.e. overdense regions) first at small scales and then at ever 
larger scales. Let us consider first more closely the temporal 
character of this evolution. It is canonical and instructive in 
the study of cold dark matter clustering in the universe to 
consider evolution from initial conditions characterized by 
power spectrum (PS) P{k) ~ (with n a constant), and 
velocities prescribed by the growing mode of linear pertur- 
bation theory. Although the "real" cosmological power spec- 
trum (in currently favoured cold dark matter models) is not 
exactly power-law it may be well approximated as such, with 
a slowly varying exponent, between n = — 1 and n = — 3 over 
the relevant range (see, e.g. lPeebled (|l993l )V The simulation 
of such "scale-free" initial conditions is expected to give rise 
to simple behaviours (which one might hope to understand 
analytically). In particular, for — d < n < 4 ( where d is the 
spatial dimension) one expects theoretically (|Peebleslll98Cll ) 
a process of hierachical structure formation which 

• is driven by linear amplification of the initial fluctua- 
tions. For power law initial conditions and the models we 
are considering, it is simple to show that this corresponds 
to the scale of non-linearity (at which perturbations are of 
order one) which grows in proportion to 



ii.(i) = exp('^^r 
' I + n 



(5) 



1 is the root of the charac- 



where A+(r) = 

teristic equation of Eq. ((3| associated to the growing mode 
of linear theory; 

• tends asymptotically to "self-similar" behaviour in the 
whole non- linear regime, i.e., the temporal evolution of all 
correlation properties can be obtained by a rescaling of 
the spatial coordinates, in proportion to the same function 
Rs{t). 



Extending analy ses in lYano fc Goudal (|l998l ) and in 
iGabrieUi et all (|2009D we have studied a range of values of 
initial PS (n = 0,2,4) and models (T = 0,1/^6,1/^2). 
The upper panel in Fig. [1] shows the measured reduced two 
point correlation function ^{x,t) as a function of x/Rs{t), 
while the lower panel shows the (dimensionless) A — kP{k) 
as a function of fc7?s(t), for the case of an initial PS with 
n = 2 and a model with P = l/y/E. The superposition of 
the rescaled PS at small k in the lower panel shows the valid- 
ity of linear theory, while the superposition in the same plot 
at larger k, and of ^{x) in the upper panel over a wide range, 
shows the development of the corresponding self-similarity 
in the non-linear regime. We observed the same behaviours 
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Figure 1. The two point correlation function (upper panel) and 
power spectrum (lower panel) as a function of rescaled spatial 
variables, for the F = l/VE model and initial conditions with 
n = 2 (i.e. EdS cosmology). The "time" variable t = log a (where 
a = 1 initially). L is the (periodic) box size. 



in all cases, except for n = 4. In the latter case it turns 
out that, as expected, linear theory breaks down at small k, 
but nevertheless self-similarity is observed in the non-linear 
regime. 

In three dimensions such self-similar behaviour has been 
observed numerically to apply in all cases which have been 
simulated. Specifically in EdS universes it h as been shown to 
apply i n the range —3 < n ^ 1 (see, e.g.. lEfstathiou et al.l 
(Il988h: IColombi et al.l (| 19961 ): I Jain fc BertschingeJ (|l998h : 
ISmith et al.l ( 20031 )). while in the static universe limit it has 
been shown to apply for n — and n = 2 (jBaertschiger et aP 
[2007, 2008i). As in all these 3D good power-law fit 

can be obtained (as in Fig. [T} to the correlation function, 



C(r) 



and to the power spectrum ( k P{k) 



the range of non-linear self-similar clustering. The exponents 
7 we measure in our ID simulations, which are reported be- 
low in further detail, show the same qualitative dependence 
on n as in 3D: in the expandi ng case 7 increases as n does so 
(see, e.g. JSmith et al.l (|2003D ). while in the static case one 
observes a mu ch smaller exponent wh i ch do es not depend 
sensibly on n (jBaertschiger et al.|[2007l . I2OO8I ). 
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4 EXPONENTS OF SCALE INVARIANT 
CLUSTERING 

Compared to 3D the power law behaviours we observe in 
the ID models extend over a much larger range of scales, 
simply because of the much greater dynamical range acces- 
sible. In Fig. [Jl for example, such behaviour extends over 
roughly four decade s in scale, compare d to a single decade 
in 3D studies (e.g. ISmith et al.l (|2003l ')'l. Such spatial res- 
olution makes it possible to robustly determine whether 
such behaviour is indicative of scale invariant (i.e. fractal- 
type) properties in the distributions, using tools such as box- 
counting methods to measure the spectrum of "generalized 
dimensions" (see, e.g. iGabrielli et al.l (|2005h ). In a recent 
study of these ID m odels with initia l cond itions like those 
we are studying here. I Miller fc Rouetl (|2010l ) have used such 
methods to analyze this clustering and have shown that, like 
in their previous studies o f slightly different initial condi- 
tions in lMiller et al.l (|2007l ). they indeed find clear evidence 
for scaling behaviours indicative of fractal clustering in this 
regime. We have performed much of the same analysis on 
our simula tions, and find results very consistent with those 
reported in iMiller fc Ro"ue^(|201(]| ): robustly measured "gen- 
eralized dimensions" in the range of scales where the power 
law behaviour in the correlation function indicates scale in- 
variance, and in particular a value of the "correlation di- 
mension" D2 « 1 — 7 as expected. 

We show now that the exponent 7 measured from ^(a::, t) 
can be very well accounted for, in the expanding models, 
using both the self-similarity of the evolution in time, and 
directly observable "stable clustering" behaviour in the sys- 
tem at sufficiently small scales. To do so let us consider the 
temporal evolution of the lower cut-off, Xmin say, at which 
the self-similar scaling breaks down. Just as can be seen in 
Fig. [T]we observe in all our simulations that (,{x,t) breaks 
from the power-law at this scale and become approximately 
flat, reaching a maximum value S,max- (In 3D simulations the 
effects of force smoothing make it impossible to follow in this 
way the development of self-similarity at small scales). We 
observe to a very good approximation, in our simulations 
with r / 0, that 



exp(— er) 



(6) 



where e is a constant which depends only on V (and not on 
n). Further we observe that the upper cut-off to power-law 
behaviour in the correlation function, Xmax, occurs (as in 
three dimensions) at a fixed amplitude in ^ slightly larger 
than unity, i.e., roughly at the non-linearity scale [defined 
by (,{xnl) = 1]. We therefore have, from self-similarity, that 
Xmax oc Rs{t). It foUows that the exponent 7, at sufficiently 
long times, should be given by 



7 = 



^ + T^A+(r) 



(7) 



To determine the exponent fully we need only deter- 
mine how e in ([111 depends on F. The relation ^ indi- 
cates the relevance of "stable clustering" : if the correlation 
up to this scale Xmin arises from structures which do not 
evolve (macroscopically) when considered in spatial coor- 
dinates rescaled in proportion to Xmin, one obtains such a 
behaviour. In three dimensions such a behaviour has b een 
proposed long ago (|Peebledll973 : iDavis fc Peebleslll977h as 
possibly valid in the strongly non-linear regime, and used 



to obtain a prediction for the exponent 7 in an EdS uni- 
verse. It amounts to supposing that if non-linear structures 
behave essentially as isolated systems (i.e. if the tidal forces 
exerted by all mass external to them are negligible), they 
will be expected to virialize and remain stable in physical 
coordinates. In the ID models we are considering the equa- 
tions of motion have not been derived starting from physical 
coordinates, but such "stable clustering" nevertheless has a 
clear meaning. Indeed in the ID model the meaning of "iso- 
lation" of a subsystem may be given more exactly than in 
3D: if particles of a given subsystem do no cross particles 
outside the system, the corresponding system indeed evolves 
independently from the rest of the "universe" (i.e. tidal grav- 
itational forces vanish in ID). To see what evolution is then 
predicted, it is appropriate to rewrite the equations of mo- 
tion for such a finite isolated subsystem using the particle 
labelling in which particles pass through one another. It is 
straightforward to show that the equation of motion for par- 
ticles belonging to such a subsystem may then be written 
(in the units we have adopted in which g = l/2no) 



d Xj ^ ^ dxj 
dr^ dr 



-g{N>{{)-N<{i))+Xi 



(8) 



where Xi is the position of the particle i with respect to the 
centre of mass of the system, and A^^ (i) (A'^^ (i)) is the num- 
ber of particles (in the subsystem) on the right (left) of the 
particle i. Unsurprisingly this is just the equation of motion 
for a finite ID self-gravitating system in an infinite space, 
with an additional damping term arising from the expansion 
and a term from the subtracted mean density. We can de- 
fine the conserved energy E associated with the two terms 
on the right hand side, and have then 

^ = -2rK (9) 
where K is the total kinetic energy, and E = K + Ugrav + 

Ubg with Ugrav = ^Y^i,jg\Xi -Xj\ aud Ubg = Y^iX'^/2. 

Now, if such an isolated subsystem is significantly overdense 
(i.e. its mean density ris is much greater than no), the time 
scale associated with the dynamics of the first term on the 
right hand side Tgrav ~ \/no/ns is much shorter than that 
associated with the other two terms, and the term Ubg is 
negligible compared to Ugrav ■ We expect this to allow the use 
of an adiabatic approximation to in which E and K axe 
replaced by their values averaged on this shorter time scale, 
(E) and {K}. On these time scales, however, we expect the 
corresponding virial relation to apply, i.e., {2K) — (Ugrav)- 
It follows that (E) = 3{K), and then, using Eq. (O, that 



(K) oc exp 



2rr 



(JJ grav^ ■ 



(10) 



Given that Ugrav is linearly proportional to the size of the 
subsystem, the size of such a virialized subsystem will be 
expected to scale in the same way. Thus, if the behaviour 
in ((6]) is indeed due to such stable virialized structures, e = 
2r/3. Substituting this relation in Q gives 



7 ■ 



2r(n-M) 



F(2n - 1) + 3VW+4: 



(11) 



^ We note that lYano &: Goudal l ll998t) give an incorrect general- 
ization to ID of the stable clustering hypothesis, assuming that 
stability will be attained in the 3D physical coordinates. 
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intial PS 


r = (thy) 


r 


= (sim) 


r = 1/V6 (thy) 


r = 


l/\/6 (sim) 


r = 1/^2 (thy) 


r = 


1/^2 (sim) 


n = 


7 = 


7 = 


0.18 ± 0.03 


7 = 1/7 


7 = 


0.14 ± 0.02 


7 = 1/4 


7 = 


0.25 ± 0.02 


n = 2 


7 = 


7 = 


0.18 ± 0.03 


7 = 1/3 


7 = 


0.35 ± 0.02 


7 = 1/2 


7 = 


0.50 ± 0.02 


n = 4 


7 = 


7 = 


0.15 ± 0.05 


7 = 5/11 


7 = 


0.43 ±0.01 


7 = 5/8 


7 = 


0.62 ±0.01 



Table 1. Values of the exponent 7(n, F) characterizing the power-law range in !;{x), as measured in simulations ("sim") and as predicted 
theoretically ("thy") by JTTJ, which gives 7 = (n + l)/(n + 7) for T = I/VG, and 7 = (n + l)/(n + 4) for T = 1/v^. 



We note that, using the above arguments, a straightforward 
generalization of the usual 3D stable clustering prediction 
(for an EdS cosmology) to the class of cosmological back- 
grounds corresponding to any (constant) F may be obtained. 
The resuh ifl 

6r(n + 3) 



73 



r(2n + 5) + ^/^2T4 



(12) 



which indeed coincides with the usual EdS prediction, 7 = 
3(3 + n)/(5 + n), for F = l/x/6. 

Shown in Table[T]is a comparison of the prediction 
with the corresponding exponents measured in our simula- 
tions. The agreement is excellent for the expanding cases, 
and indeed in these cases the value e measured directly 
agrees well with the predicted e = 2r/3 0. In the static 
case the prediction leads to the badly defined result 7 = 0, 
while we observe a small exponent 70 ~ 0.16 which does 
not appear to depend on n. We observe in this case that, 
rather than ((S)), we have an Xmin which remains approxi- 
mately constant (i.e. the expected F = limit), but a ^max 
which grows. Clearly the stable clustering hypothesis is in- 
deed only an approximation, to which corrections (due to 
the fact that structures are not, of course, strictly isolated) 
become most important, as one would expect, when we go 
to the limit F — >■ in which virialized structures no longer 
decrease in size. 



5 HALOS AND VIRIALIZATION IN THE 
SCALE INVARIANT REGIME 

Let us now consider the observed clustering in terms of "ha- 
los", in analogy to how this is done in cosmological sim- 
ulations. We thus use a friend-of-friend algorithm, with a 
variable linking length A (which we define in units of the 
initial lattice spacing £). In ID this is equivalent to sepa- 
rating the system into structures separated by voids of size 
greater than or equal to A. In the range of scale-invariant 
non-linear clustering, i.e., for Xminjl < A < xnl/1-, we ob- 
serve that the algorithm selects out structures with a broad 
range of sizes from our simulations, peaked around a value 
which scales with A. However, the fact that the distribution 
of sizes is quite large, while the range of scale invariance 



^ Virialization in three dimensions leads to {E) = —{K), so that 
{Ugrav) oc exp(2rT). It foUows that x^i^ oc cxp(— 2rT), and 
then, with ^max c< x^i„ and the 3D expression for Rs(t), the 
result follows. 

For the case F = l/\/6 our me a sured exponents agree with 
those reported by lYano &: Goudal lll998l ). Our stable cluster- 
ing result e = 2F /3 explains also the exponent measured by 
lAurell et al.l ||2003| ) directly for an isolated structure. 



in any simulation is of course quite limited, make it difficult 
(even in ID) to use such halos as a tool for characterizing this 
regime. We thus analyse here the case of an initial PS with 
n = 4 in which this range is greatest — a little more than 
four orders of magnitude in scale between Xminji ~ 10^^"'' 
and XNh/t ~ 10^'^ in the most evolved configuration. 

Shown in Fig. [2] is the measured distribution of virial 
ratio as a function of A in this simulation, i.e. the distribu- 
tion of the values of the ratio 2K /Ugrav measured directly 
in each halo, taking the velocities with respect to the centre 
of mass of the corresponding halo. We observe an approxi- 
mately stable distribution peaked about unity in the range of 
A where the structures selected out fall mostly in the length 
scales of scale invariance, while at the smallest scales no such 
tendency is observed. On the contrary the virial ratio tends 
to be very large, which can be interpreted as due to the fact 
that the algorithm is selecting structures inside the scale 
Xmin- if the distribution at this scale is constituted of ap- 
proximately smooth virialized structures of size Xmin, parts 
of these structures will have a high virial ratio (because K 
scales more slowly with size than Ugrav in a smooth struc- 
ture). 

We have studied quantitatively the stability of these 
measured distributions in the range of scale invariant clus- 
tering using a nonparametric Kolmogorov-Smirnov (K-S) 
test. The K-S statistic quantifies a distance between the em- 
pirical distribution functions of the different samples, and 
the null distribution is calculated under the null hypothe- 
sis that the samples are drawn from the same distribution. 
Calculating the corresponding p- values of this test, we find 
that the null hypothesis is indeed not rejected in the range 
of scale-invariance, while, on the other hand, it is clearly 
rejected outside this range. 

Our interpretation of these results is that one can effec- 
tively decompose the distribution into a collection of struc- 
tures which are, approximately, virialized. This is true pro- 
vided these structures are defined at any scale in the range 
of scale-invariant clustering. 



6 DISCUSSION 

In three dimensions the stable clustering hypothesis has 
been used widely as a reference point, and indeed much used 
phenomenological models such as the formalism proposed by 
iPeacock fc Doddsl (|l996l ) incorporate it into the modeling 
of the non-linear power spectrum obtained from numerical 
simulations. However, in parallel, results of major simula- 
tions o f power law i nitial conditions (e.g. lEfstathiou et al.] 
(|l988l ): ISmith et al] (|2003l )) led to the conclusion that, al- 
though the measured exponents in the correlation function 
showed a behaviour roughly consistent with its predictions 
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Figure 2. Normalized distribution of the virial ratios of halos selected with the indicated values of A (linking length in units of lattice 
spacing), for the case n = 4 and F = l/\/6 at t = Ina = 22. 



(albeit with some significant deviations (|Smith et aLlbOO^ )). 
there was no evidence for the "clustering hierarchy" which 
Peebles had a rg ued would be ass o ciated with it (see, e.g. 
IPeeblesI (|l980l ')'). lEfstathiou et al.1 (|l988l 'l. notably, explic- 
itly excluded such a "fractal" description and found evi- 
dence instead for the validity of a description in terms of 
"smooth non-linear clumps". These latter are the precur- 
sors of the "halos" of halo models, which have become the 
standard phenomenological description of the matter distri- 
bution in cold da r k mat ter cosmologies (for a review, see, e.g. 
ICoorav fc ShethI (|2002h l. The matter distribution is then ap- 
proximated as a collection of isolated (and thus virialized) 
spherical structures with smooth density profiles. Further, 
on the basis o f exten sive numerical study following that of 
iNavarro et al.) Ewh . the latter are widely believed to be 
characterized well by "universal" exponents (i.e. indepen- 
dent of initial conditions and cosmology) . 

In one dimension the strongly non-linear regime is 
truly scale-invariant (for power law initial conditions) in a 
range which grows monotonically in time. The associated 
distribution of matter is intrinsically lumpy or grainy down 
to the lower cut-off scale Xmin- indeed the very meaning 
of such scale invanance is that there are no characteristic 
scales available to define smoothness. To illustrate this we 
show in Fig. [3] the spatial density at various levels of de- 
tail in a typical evolved configuration. Such distributions 
are most naturally described with the instruments of frac- 



tal (or, more generally, multi-fractal) analysis developed for 
this purpose in condensed matter and statistical physics (see 
e.g. iGabrieUi et"al] (120051 )1. As we have just illustrated in 
the previous section, halo type descriptions may, of course, 
also be employed to describe them. The halos so defined are, 
however, very different to those described by 3D halo mod- 
els: these ID halos are not smooth. Further have they no 
intrinsic size themselves, but are defined only with reference 
to an arbitrary chosen scale. The study of the virial ratios 
we have presented indicates, however, that such halos can 
be considered as entities with a dynamical relevance, as they 
show a clear tendency to have a virial ratio of order unity 
(which is the behaviour of an isolated structure). The clus- 
tering in the non-linear regime can thus be considered as a 
concrete realization of the qualitative picture of a "cluster- 
ing hi erarchy" originally envisaged by Peebles (e.g. [Peebles! 
(Il980l )) as resulting from stable clustering. The stable clus- 
tering hypothesis we have described above, however, is actu- 
ally subtly different from the original one: we assumed only 
that stable clustering applies below the scale Xmin marking 
the lower cut-off to the scale invariance, and not necessarily 
to the strongly non-linear regime as a whole. Thus we as- 
sumed only that stable clustering applies at an ultraviolet 
scale fixed by the resolution of the simulation (or, physically, 
by the scale at which the very first structures form). The 
"statistical virialization" we have observed using the halo 
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Figure 3. Density field obtained at t = logo = 14 starting from an initial condition with n = 2, in the "EdS" model with T = 1/^6. 
The first panel shows the whole box (i.e. length units are such that L = N = lO''), while each subsequent panel shows a spatial "zoom" 
on a region of size one tenth that shown in the previous panel. The resolution has been increased (i.e. the bin size has been decreased) 
in the last three plots in order to reveal the clumpy nature of the distribution at these scales 



analysis, on the other hand, apphes at scales above Xmin 
and across the range of the scale invariant clustering. 

There are clearly two possible conclusions one can draw 
from this analysis: 

• A. These ID models produce non-linear clustering 
which is qualitatively different in its nature to that in 3D, 
or 

• B. The spatial resolution in 3D simulations up to now 
has been too limited to reveal the nature of clustering in 
cold dark matter cosmologies, which is correctly reflected 
(qualitatively) in the ID simulations. 

We believe that, despite the impressive computational 
size and sophistication of 3D cosmological simulations, con- 
clusion B may well be the correct one. The very largest 
modern studies in a cosmological volume access roughly two 
decades in scale in the non-linear regim^B while reference 
studies in the literature of power law initial conditions in 



^ F or example, the "millennium" simulation llSprineel et al.l 
l2005h of the KCDM cosmology, has, at redshift zero, ^ = 1 at 
about 2h~^Mpc, an initial mean interparticle I ^ 250h~^ kpc 
and a force smoothing e = 5h~^ kpc. It is usually supposed that 
it is the latter which sets the lower bound on resol ution, but we 
note that it has been shown by several studies fsee lMelott et ahl 



EdS cosmology jEfstathiou et al.l 1 19881 : ISmith et all |2003| ) 
measure the crucial power-law behaviour in the correlation 
function (or the PS) over at most one decade. If we were 
to perform our ID simulations at c omparable resolutio n to 
large cosmological simulations like ISmith et al.l (|2003l ). we 
would certainly have great difficulty in establishing the scale- 
invariant nature of the strongly non-linear clustering arising 
from power law initial conditions. Although halos defined 
exactly as in three dimensions might look clumpy, an ap- 
proximately smooth profile could be determined for them 
if they were averaged (as they can be in three dimensions 
when spherical symmetry is assumed). Higher resolution 
3D simulations of smaller regions have shown over the last 
decade that there is in fact much more substructure ("sub- 
halos") insid e halos t han was originall y antici p ated (see, e.g., 
iMoore et all (|l999l ): iDiemand et~all (|2005l ): iGoerdt et al] 
(20o3)), and have even more recently described several lev- 



els of such sub s truct ur e ("subhalos o f subh al os", see e.g. 



iDiemand" et all (l2008l); ISpringel erall (|2008D : IStadel et al.l 
( 20091 )). [pTemand et al] (|2008l ) everTuse the term "fractal" 
to describe (qualitatively) the real space structures, while 



lll997)) : lRomeo et al.l l f2008l) : |jovce et al.l | |2009|) ) that, at the very 
least, precision may be compromised below the scale I. 
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IZemp et all (|2009l ) describe the structure of halos in phase 
space as "intrinsically gr a iny". We not e that other authors 
(see, e.g. JValagea^ (|l999l ): lGait j (j2007l )l have previously ar- 
gued for similar conclusions on the basis of analyses of 3D 
simulations. 

Let us consider nevertheless one possible consideration 
in favour of (the more conservative) conclusion A. In the 
expanding (i.e. damped) ID models, the stable clustering 
prediction (|lip fits the measured exponent s extremely well. 
Early 3D studies for EdS cosmologies (e.g. lEfstathiou et al.] 
([1983)) measured exponents roughly consistent with the sta- 
ble clu stering prediction, but later studies (e.g. ISmith et al.l 
have found significant disagreement. This disagree- 
ment is attributed to physical mechanisms which cause the 
fundamental assumption of stability to be violated — by 
the evident fact that there are interactions between "halos" , 
which can even lead to their merging into single structures. 
We have noted that in one dimension tidal forces vanish, and 
structures can interact only when they actually physically 
cross one another. While merging may occur, it may be that 
it is a less efficient process than in three dimensions. Thus 
the excellent agreement in the ID models compared to EdS 
may perhaps be attributed to the fact that these models 
probably represent poorly the role of such physical effects. 
The essential question, however, is not whether these effects 
play a role and can lead to deviations from stable clustering, 
but whether such effects can qualitatively change the nature 
of clustering, destroying scale invariance by smoothing out 
the distribution on a scale related to the upper cut-off to 
scale invariance. Our study of the case F = suggests that 
the answer is negative. The prediction of stable clustering 
does not work in this case, and like in three dimensions, one 
obtains a small value of the exponent which does not sensi- 
bly depend on n. The physical reasons why the exponent is 
close to, but different to, the stable clustering prediction are 
a priori the ones just cited. The analysis of Miller et al. of 
this case, which we have rechecked and confirm, finds nev- 
ertheless that the distribution is scale invariant. Further, as 
we have mentioned, the lower cut-off Xmin remains constant 
as in the stable clustering hypothesis, of order the initial 
lattice spacing (and unrelated to the upper cut-off). 

These results on ID models suggest directions for 3D in- 
vestigations which might establish definitively the correct- 
ness of conclusion B. We note, for example, that the ID 
models lead one to expect that the exponents derived phe- 
nomenologically to characterize the highly non-linear den- 
sity field inside smoothed halos (i.e. the "inner slope" of ha- 
los) should be closely related to the exponent 7 determined 
from the correlation function. Indeed — in the approxima- 
tion of a simple fractal behavior in the strongly non-linear 
regime, which the spectrum of m ulti-fractal exponents mea- 
sured in iMiller fc Rouetl (|201Cll l suggests should be quite 
good — the mean density about the centre of such halos 
will decrease just as about any random point, i.e., with the 
same exponent 7. Despite the contradiction with the widely 
claimed "universality" of such exponents in halos profiles, 
such a hypothesis cannot currently be ruled out, as the de- 
termination of such exponents is beset by numerical difficul- 
ties (arising again from the limited resolution of numerical 
simulations). In a stu dy of halo profiles obtain ed from power 
law initial conditions iKnoUmann et aL I (|2008D show explic- 
itly that the results for the halo exponents depend greatly 



on what numerical fitting procedure is adopted. While one 
procedure gives "universality" (i.e. exponents independent 
of n) , a different one favors clearly steepening inner profiles 
for larger n. Indeed we note that the numerica l valu es for 
the inner slopes obtained by iKnoUmann et all (|2008l ) are, 
for the larger n investigated, in quite good agreement with 
the exponent predicted by stable clustering. 

Our results here are for power law initial conditions, 
but we could equally use the model to study both initial 
conditions and an appropriately modified time-dependent 
damping rate mimicking the ACDM model in one dimen- 
sion. The exact scale invariance in the strongly non-linear 
regime would certainly be broken in this case, and the var- 
ious characteristic scales introduced will be imprinted on 
the clustering. There would thus not be a single correlation 
exponent, but a slowly varying one. One would certainly ex- 
pect, however, the qualitative nature of the clustering to be 
unchanged, just as in three dimensions clustering appears 
to be qualitatively the same in "scale-free" cold dark matter 
cosmologies and ACDM. Further, our considerations here 
are strictly relevant only to dissipationless cold dark matter 
simulations. If the initial conditions are "warm" or "hot", 
or if other non-gravitational interactions are turned on, the 
associated physical effects will tend to smooth out the mat- 
ter distribution up to some scale (and thus destroy the scale 
invariance up to this scale). Nevertheless, if the conclusion 
B is correct even for this idealized case, it is likely to have 
very important observational implications relevant to test- 
ing standard cosmological models. Intrinsically clumpy or 
grainy halos lead, for example, to ver y different prediction s 
for dark matter annih i lation (see, e.g. iGoerdt et"aL 
iDiemand et~aLl (|2008h : lAfshordi et al.l (|2010l ')')'). Concerning 
the compatibility of such a distribution of cold dark mat- 
ter with observations of the distribution of visible matter at 
sub-galactic scales — which would be expected to lead to 
an amplified "missi ng satellite" problem — we n ote that re- 
cent studies such as lWadepuhl fc Springell (|201Cf ) show that 
solid conclusions in this respect will be reached only on the 
basis of a much improved understanding of the processes in- 
volved in star formation. At larger scales the possible link to 
the striking power-law behavior which characterize s galaxy 
correl a tions over several d e cades i n scale (see, e . g.. iPeebled 
(|l974h : [Svlos Labini et al.l (|l998l ): iMasiedi et al.l (|2006l 'l') — 
which was the motivation for original work using stable clus- 
tering t o explain such behaviour from power law initial con- 
ditions (|Peeblesl [l97i ) and is naturally interpreted as in- 
dicative of underl ying scale invariance in the mat t er dis - 
tribution (s e e, e.g . ISvlos Labini et all (|l998l ): ICaitd (|2007h : 
lAntal et al.l (|2009l )') — is intriguing, and will be discussed 
elsewhere. 

We acknowledge useful discussions with Andrea 
Gabrielli, Jean-Michel Levy, Francesco Sylos Labini and 
Tirawut Worrakitpoonpon. We thank the anonymous ref- 
eree for several very useful suggestions. 
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